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We consider quantized vortices in two-component Bose-Einstein condensates and three-component 
Fermi gases with attractive interactions. In these systems, the vortex core can be either empty 
(normal in the fermion case) or filled with another superfluid. We determine critical values of 
the parameters - chemical potentials, scattering lengths and, for Fermi gases, temperature - at 
which a phase transition between the two types of vortices occurs. Population imbalance can lead 
to superfluid core (coreless) vorticity in multicomponent superfluids which otherwise support only 
usual vortices. For multicomponent Fermi gases, we construct the phase diagram including regions 
of coreless vorticity. We extend our results to trapped bosons and fermions using an appropriate 
local approximation, which goes beyond the usual Thomas-Fermi approximation for trapped bosons. 



; I. INTRODUCTION 

Properties of quantized vortex cores in superfluids have 
been actively researched for many years. For example, it 
was realized [l| that in type II superconductors there 
are low-energy states bound to the core and their ef- 
feet on the local density of states was studied experimen- 
tally [H and theoretically 3 . More recently, observations 
of quantized vortices provided key evidence for super- 
fluidity in both single component Bose-Einstein conden- 
sates (BECs) and two-component superfluid Fermi gases, 
where vortex cores are detected as regions of suppressed 
particle density [1, JH|- The core states in Fermi gases 

[ were considered in [6|. 

The situation in unconventional superfluids is more 
complex. For instance, in superfluid ^He-B the vor- 
tex core is in a ferromagnetic superfluid state fj\. In 
high-temperature superconductors with d-wave order pa- 
rameter an s-wave component must be present in the 
core @. In a color superconductor - high-density, low- 
temperature quark matter - there are many distinct 
fermion species (components) that can pair up, leading to 
two types of vortices - Abelian and non- Abehan [3]. sim- 
ilarly, spinor atomic BECs having several bosonic com- 
ponents host, in addition to vortices, other types of topo- 
logical excitations such as hedgehogs and skyrmions [13] ■ 
This cannot happen, however, if the symmetry between 
the components is explicitly broken, in which case only 
vortices are possible. 

Consider, for example, a three-component Fermi gas , 

' which can support three distinct superfluid states [111 . 
IT 2] , with a vortex in superfluid Si . The core of the vor- 
tex can be in the normal state, as for a two-component 
Fermi gas, or it can condense into one of the two other 
superfluid states - see Fig. [TJ In view of current efforts 
to achieve superfluidity in three-component Fermi gases 
[T3I [l^ . it is important to understand which of these 
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FIG. 1: Order parameters for normal-core [panels (a) and (c)] 
and superfluid-core [(b) and (d)] vortices in multicomponent 
BECs and Fermi gases. Colors distinguish different superflu- 
ids. Top row: qualitative behavior of the order parameters 
in a plane of a trapped superfluid. Color intensity is pro- 
portional to the local value of the order parameters. Bottom 
row: profiles of the condensate wave functions in the absence 
of trapping. 

scenarios is realized and/or how to drive the transition 
between normal- and superfluid-core vortices. Similarly, 
competition between empty- and filled-core vortices takes 
place in a two-component BEC (2BEC). In fact, fiUed- 
core vortices have been experimentally observed fl6j , 
and their properties are in the focus of numerous the- 
oretical studies [13, in, [S S mi - see Ref. M for a 
review. 

In this paper, we study vortices in 2BECs and three- 
component Fermi gases. We show that both empty- 
(normal-) and filled-(superfluid-)core vortices can be re- 
alized depending on the parameters characterizing the 
bosonic (fermionic) system: chemical potentials, scatter- 
ing lengths, and temperature in the case of fermions. We 
first focus on 2BEC at zero temperature and derive, us- 
ing Gross-Pitaevskii equations, the critical relationship 
between the parameters of the system which determines 
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whether the vortex is empty or filled. Fermions at suf- 
ficiently high temperature can be analyzed in a simi- 
lar manner with the help of the Ginzburg-Landau ex- 
pansion. This enables us to find a condition for the 
normal- to superfluid-core transition for vortices in a 
three-component Fermi gas. In particular, for balanced 
systems with equal populations of three fermion species, 
we explicitly obtain the critical temperature for this tran- 
sition. We first consider bosons and fermions in the ab- 
sence of an external potential, then extend our results to 
trapped systems. We further discuss the validity of the 
widely used Thomas-Fermi and local density approxima- 
tions for analyzing vortex cores in trapped condensates. 




II. TWO-COMPONENT BOSE SYSTEM 

The low temperature properties of a BEG are well de- 
scribed in terms of a condensate wave function 5* which 
obeys the Gross-Pitaevskii equation, which captures the 
effects of interactions at the mean field level. First, let 
us treat a 2BEG free of the trapping potential. The cou- 
pled Gross-Pitaevskii equations for the wave functions 
^i, i = 1,2, of the two condensates in a stationary state 
can be obtained by minimizing the following energy den- 
sity functional [To ]: 
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where rm and fii are the atomic masses and chemical 
potentials, respectively. The coupling constants Uij are 
related to the (positive) scattering lengths 



Uij — 2Traij {l/rrii + \/m-j) . 



(2) 



We consider the case in which the two BECs phase 
separate (i.e., — W11W22 > 0) and assume that there is 
a vortex line along the z axis in condensate 1. Then the 
wave function of the condensate has the form 



*i(r,' 



(3) 



where ^1 — l/^/2mi/ii is the healing length and r = 
(r, z). The profile function f{x) is the solution of the non- 
linear differential equation (primes denote derivatives) 



f f 



(4) 



with the boundary conditions f{0) — and f{x) — > 1 as 
X 00. 

Condensate 2 fills the core of the vortex in condensate 
1 when this is energetically favorable. At the transition 
from and empty (4*2 = in the core) to a filled core 



FIG. 2: Ground state energy eo of Schrodinger Eq. (O as a 
function of the inverse mass /3,.. Solid line: numerical calcula- 
tion. Dotted (dashed) line: analytical approximation Eq. ([8]) 
for /3r ^ 3 (/3r ^ 3). Inset: relative deviations between the 
approximations and the numerical results. 



(^'2 7^ 0), 5*2 is infinitcsimally small. The energy differ- 
ence between filled and empty core states to the lowest 
order in ^'2 is 



SE 



1 



2m2 



IV^'2 



Ul2|*in«'2p -Ai2|^'2p 



(5) 

where is the size of the system in the z direction. We 
can neglect the effect of vj/j on 'Si, since it is of higher 
order in ^'2. This equation can be rewritten as 



SE = L,ei 
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Mil 



eo(/3r)/^i - ^J■2 



(6) 



where eo is the (dimensionless) ground state energy of 
the following two-dimensional Schrodinger equation: 

- PrV^^2ip, 0) + /'(P)*2(P, 0) = eo*2(p, 0) (7) 

with p = r/^i and the "inverse mass" being Pr — 
miuii/m2Ui2- The analysis reported in Appendix |A] 
shows that a good approximation for the monotonically 
increasing function eo(/3r) is given by (see also Fig. [21) 



eoiPr) 
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/3r <3, 
>3, 

(8) 

where 7^ is Euler's constant and cq = /'(O) ~ 0.58319. 

As an example, consider a system with mi = m2 and 
Ml ^ M2- Then Eqs. ([6]) and ([8|) show that for a22 < 
012 < flu a coreless vortex in supcrfluid 1 is more stable 
than a coreless vortex in superfluid 2, since the latter 
state has higher energy. This is in agreement with the 
experimental findings of Ref. [l^. 

More generally, the sign of the term in square brackets 
in Eq. ^ determines the stable core state. If it is pos- 
itive, it costs energy to fill the core. Hence the core is 
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empty if 
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In other words, Eq. ([9]) defines a surface in the space 
of parameters that separates the regions of empty and 
filled core vorticities. For example, we see that popula- 
tion imbalance, which modifies the ratio /ii//i2, impedes 
coreless vorticity in one of the condensates while favoring 
it in the other, as seen by exchanging 1 <-> 2. Analogous 
conclusions hold for the ratio between intra- and inter- 
species scattering length and the ratio of masses. Inter- 
estingly, while our calculations are performed in the ther- 
modynamic limit, qualitatively similar conclusions were 
reached in Refs. llTHlsL where stability conditions were 
derived for vortices in a small trapped condensate. 

Now, let us analyze the effects due to an external po- 
tential. For simplicity, we consider a spherically sym- 
metric harmonic trap V{v) and assume, without loss of 
generality, that condensate 1 occupies its center, while 
superfluid 2 forms a shell around it. Vortices in trapped 
2BEC have been previously studied in Refs. [l^, [lo, [U 
using the Thomas-Fermi approximation (TFA). In this 
approximation only the angular part of the kinetic en- 
ergy terms in Eqs. ([1]) and (O is kept. Then, the wave 
functions ^'i, i = 1,2 of the condensates are obtained 
from those in the absence of the trapping potential via 
the replacement /i^ = Mi(0) — ^(r), i.e., the 

trapped 2BEC is taken to be locally uniform. This is 
justified in large condensates, Ri 3> ^i, where Ri and 
are condensate sizes and healing lengths, respectively. 

Neglecting the radial parts of the kinetic energy, how- 
ever, is a good approximation only when ^'^ varies over 
distances much larger than ^i. As we have seen above, 
this is not the case close to the vortex line, where both 
^1 and \E'2 vary on the scale ^i, see Eqs. (|4|) and d?]). We 
therefore expect the TFA to break down in determining 
the state in the vortex core. Indeed, one of its artifacts 
is that is identically zero in a finite region inside the 
core [l^, im . This implies that the second term in Eq. ^ 
vanishes too, while the last term always makes it energet- 
ically favorable for the otherwise empty core to be filled 
by the second superfluid. Thus, the question whether 
the core is empty or filled cannot be resolved within the 
TFA, prompting arbitrary assumptions of filled [l^ [l^l 
and empty [21| core vortices in the literature. 

This question, as we now show, can be accurately 
answered for large condensates by combining our ap- 
proach with the local uniformity assumption discussed 
above. For example, let us work out the condition under 
which an empty core vortex is realized in superfluid 1 for 
TOi = 7712. The separation of scales, Ri ^ ^i, allows us to 
describe the vortex profile by Eq. ([3|) with fii — > /x^ (r) . In 
other words, we neglect gradients of ^i{r) but not of f{r). 
This approximation is valid at distances L ^ from the 
interfaces between the condensates [l^]. In a plane per- 
pendicular to the vortex line, the wave functions^'i looks 
like Figs.[ljc) and (d), except for the curvature imposed 



by the trapping potential. vE'i(r, (j>) changes rapidly near 
the core as r varies over distances r ~ and is smooth 
on this scale when moving along the z axis (vortex line). 

Therefore, the stability of an empty core vortex for any 
z is determined by Eq. ([9]) with /i^ — > iii{rc,z), Tc being 
the position of the core. We obtain 



fJ-2{rc,z) ^ ai2 ^ / Qii 
/ii(rc,z) an °\ai2 



(10) 



where we took into account 



77il 



= 777 and 



7772 = ™ 



A-Kaij/m - see Eq. ©. Further, one can show that 
fj.2{0) > Mi{0) is necessary for superfluid 1 to occupy 
the center of the trap [U, [2^. It follows that the 
left-hand side of Eq. (jTU]) is a monotonically increas- 
ing function of \z\, i.e., it reaches its maximum at the 
interface between the two condensates. The position 
Tin of the interface is determined from the condition 
M2(rm)/Mi(rm) = \/ 022 / Qj 1 obtained by equating the 
pressures on the two sides [111, [23| , so that Eq. (fTOjl be- 
comes 



a22 ai2 / an 

< eo 

an ail \ai2 



(11) 



If this inequality holds, it costs energy to introduce the 
second condensate everywhere along the vortex core and 
the empty vortex is the stable state. Otherwise, the core 
is partially or fully filled. Since the left-hand side of 
Eq. (|10p is larger at the interface, when varying the scat- 
tering lengths the core will be filled starting from the 
interface between the condensates towards the trap cen- 
ter. Us ing Eq. ([8]) and the experimental values reported 
in Ref. [l^, we find that condition pT|) is violated in this 
experiment. This invalidates the approach of Ref. 21 for 
the description of this experiment based on an empty 
core assumption. 

A few comments are in order about the present deriva- 
tion. We have not used anywhere the condition 0^2 — 
011^22 > that ensures phase separation in the ab- 
sence of external potentials. For trapped 2BEC, phase 
separation can occur as long as 012 > 0, although if 
a^2 — 011022 < a coexistence region is present |2Q] . 
In this case our analysis applies at distances L ^ from 
the coexistence region, but the inequality pTjl is always 
violated since eo < 1 - see Eq. ([5]). This means that 
empty vortices are possible only in phase-separated con- 
densates. More generally, as discussed above, to obtain 
Eq. (fTTjl we use locally the result derived in the ther- 
modynamic limit. This approximation is valid assuming 
slow variation of the wave functions along the vortex line. 
This assumption breaks down near the interface between 
the two phases, and even in fully phase-separated 2BEC 
the regions where the vortex core meets the interface are 
beyond the reach of the present approach. 
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III. THREE-COMPONENT FERMI SYSTEM 

The treatment presented above can be applied to the 
study of a vortex in a muhicomponent Fermi system at 
sufhciently high temperature. Consider, in particular, 
a three-component system at weak coupling near second 
order phase transition lines. The Ginzburg-Landau (GL) 
expansion for the thermodynamic potential Q is [12i] 



- / I A |2 , 



|A,|4 + iy|VA,p 



-/3l2|AlnA2|^ 



(12) 



where Jljv is the normal-state potential, v is the density 
of states at the Fermi energy, vp is the Fermi velocity, and 
Ai is the order parameter describing pairing of particles 
belonging to species j and k with j, k, and i all different. 
The coefficients ai, (3ij in Eq. p2|) are 
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0, 



where /x^ are the chemical potentials of the three species, 
ij) is the digamma function, j3ii is obtained from (3ij in the 
limit hj —f hi, and Tc^ (Tcj) is the critical temperature 
for the normal-superfluid transition for a two-component 
gas of species 2 and 3 (1 and 3). With no loss of gen- 
erality we assume Tc^ > Tc^- For simplicity, we neglect 
the weakest of the three possible interactions, say be- 
tween species 1 and 2 (recall that due to Pauli exclusion 
principle only interspecies scattering is possible in the s- 
wave channel). Then, there are only three main phases 
of the three-component system - the normal phase N 
(Ai = A2 = 0), superfluid Si (Ai ^ 0,A2 ^ 0), and 
superfluid ^2 (Ai = 0, A2 7^ 0). The applicability of the 
GL expansion reqmres T > 0.56Tc, and \h^\/4:^TT < 0.30. 
In this regime Pij > 0, indicating repulsion (this is a con- 
sequence of Pauli exclusion), and (3^2 — > 0, leading 
to phase separation between the superfluid states. The 
first-order S1-S2 transition is accessible within the GL 



description when Tc^ — Tc^ ^ Tcg and Tc^ 



see Ref. J2 for more details. 

Comparison of Eqs. ([T]) and (|12p shows that, up to a re- 
definition of the various parameters, the energy argument 
discussed in Sec.|ll]can be applied to the thermodynamics 
of the three-component Fermi gas. For example, a vortex 
in superfluid 5*1 is described by the order parameter [cf. 
Eq. ®] 



Ai(r, 



-ai 



(14) 



with the coherence length ^1 = vf\//3ii/6(— ai). 

Repeating the previous analysis for such a vortex we 
find the following condition for the second-order phase 



transition between normal and superfluid core (or equiv- 
alently between standard and corcless vortex): 



a2 



P12 (P22 



-0, 



(15) 



where the function eo(x) is given by Eq. ([8]). The first 
term in Eq. (jlSp originates from the energy gained by 
condensation of the originally uncondensed species. The 
condensation takes place in the core region, where pairs 
of the first superfluid are broken due to vorticity. The 
second term is due to the energy costs associated with 
deforming the order parameter A2 [kinetic energy in 
Eq. ([7])] and the repulsive interaction between the two 
superfluids (potential energy). 

At a given temperature T < Tc^, Eq. p5)) determines 
the area in the hi-h2 space where the core is filled. It 
can be satisfied only in the central region of this space 
where both coefficients ai^2 are negative, i.e., when con- 
densation is in principle possible in both channels. We 
show two examples of phase diagrams in Figs. [3] and H) 
At a temperature close to Tc^ (Fig. ^ the core is filled 
only in small regions around the first order phase transi- 
tion between the two superfluid states. As the tempera- 
ture is lowered (Fig. 5]) the regions where superfluidity is 
possible expand, and so do the regions of coreless vortic- 
ity. In particular, as the temperature decreases, coreless 
vortices become possible even in balanced systems with 
equal populations {hi = /12 = 0). In this case the coef- 
ficients Pij all coincide and from Eq. we find that 
the temperature To for the onset of coreless vorticity is 
related to the superfluid critical temperatures as 



7 = 



eo(l) 
l-eo(l) 



2.92. 



(16) 



Above To the core is in a normal state, while for T < To 
it is superfluid. 

Eq. (fT5|) can be used to study trapped Fermi gases 
within the local density approximation (LDA) so long 
as condensate size R is large compared to the coherence 
length ^1. Then, we can simply substitute local values 
of the critical temperatures Tc- (r) into Eq. (note 
that chemical potential differences hi are position inde- 
pendent within LDA). The condition i? ^ ^1 is satisfied 
if the total particle number A''^ is large enough. For a 
somewhat weak interaction fcF|as| — 0.5 the condition 
on the particle number is Nt ^ 10^, see Ref. [l^. We re- 
mind, however, that for imbalanced systems the presence 
of domain walls could require larger particle numbers for 
the LDA to be valid. 

The above calculations are valid in the weak coupling 
regime near second order phase transitions. The latter 
requirement is satisfied for all relevant chemical potential 
differences if Tci — Tc^ ^ Tc^ and T is close to Tc^ ■ Relax- 
ing these conditions will not alter the qualitative picture, 
although it will affect the quantitative results. For exam- 
ple, if the difference between the critical temperatures, 
— Tc2, is large, we expect the actual onset tempera- 
ture To to be smaller than that predicted by Eq. : at 
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FIG. 3; High temperature (T — 0.93rc2) phase diagram for a 
three-component Fermi gas with Tci/7c2 = 1-04 in the plane 
of chemical potential differences hi. Solid lines: second-order 
normal-superfluid phase transitions. Dashed lines: first-order 
transition between the two superfluid states. In dark grey are 
the regions where the core of a vortex in the superfluid is 
filled by the S2 superfluid. Light grey: 5*1 superfluid flUs the 
52 vortex core. 
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FIG. 4: Phase diagram as in Fig. [3] but for a lower tempera- 
ture (r = 0.85rc2). Note that now the coreless vortex state 
is present even for balanced gases with hi = 0. 



low temperature T the order parameter in the core 
region rises on a length scale of the order of the Fermi 
wave length [l^ rather than the much longer (at weak 
coupling) coherence length. At lowest order this is equiv- 
alent to an increase of cq in, e.g., Eq. ([8]), which would 
lead to a higher value for eo(l) in Eq. (|16|1 and hence a 
lower onset temperature. 



IV. SUMMARY 

We considered the empty- and filled-core vortex 
states in two-component BECs described by the Gross- 
Pitaevskii equations, as well as the normal- and 
superfluid-core vortices in three-component Fermi gases 
using a Ginzburg-Landau approach. In the absence of 
external potentials, we derived the conditions [Eq. ([9]) 
and (fT5| for Bose and Fermi systems, respectively] that 
determine the transition between the two states in terms 
of chemical potentials, scattering lengths and, for Fermi 
gases, temperature. In particular, we obtained a simple 
expression, Eq. ()16p . for the onset temperature of core- 
less vorticity in the population-balanced Fermi gas. We 
also established when the vortex core remains empty in 
trapped Bose systems, see Eq. ([TT|). using a local ap- 
proximation which goes beyond the usual Thomas-Fermi 
Approximation. We showed that, for equal masses of the 
components, empty vortices are possible only in phase- 
separated 2BECs and that in partially filled vortices the 
superfluid part of the core is in the region closer to the 
interface between the two condensates. We similarly 
applied our findings to trapped multicomponent Fermi 
gases within the Local Density Approximation and dis- 
cussed the limits of validity of this approach. The de- 
tailed study of the superfluid core is left to future work, 
as is the extension of our results to Fermi-Bose mixtures. 
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APPENDIX A: CALCULATION OF THE 
GROUND-STATE ENERGY 

In this Appendix we present the calculation of the 
ground-state energy cq of Eq. ([7]). The potential term 
f^{p) is not known explicitly, except for its behavior at 
small and large x. Indeed, near the origin a power series 
for f(p) can be found in term of one unknown parameter, 
Co, which is calculated using the boundary condition at 
infinity. The first few terms in this expansion are: 



f{p) = coP ■ 



Co(l 



192 



+ 0{p'). (Al) 
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The coefScient cg — 0.58319 can be evaluated with great 
precision either numerically or with analytical methods 
[25|. As p — !■ +00 the asymptotic expansion is: 



1 



9 



(A2) 



By appropriate rescalings, we show that this is sufRcient 
to obtain analytical estimates for eg at small and large 
(3r which, moreover, provide accurate estimates even at 
intermediate values. 



1. Small 

The limit corresponds to a particle with large 

mass; therefore, its ground state wave functions does not 
extend far from the origin. This enables us to calculate eg 
perturbatively, using the harmonic oscillator as starting 



- «l/4^-l/2^ 



Eq. © 
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point. Indeed, after the rescaling p — (3r 
becomes 

with eo — ^al\fWrCQ and Vp^{x) the perturbation poten- 
tial. The latter is defined as the potential term of Eq. ([7]), 
f^{p), minus the harmonic part. From Eq. (|A1[) we find 
the first few terms in the small-x expansion: 
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(A4) 

Consistently with the truncation of the potential Vp^, 
we calculate eq to second order in the small parameter 
via standard time-independent perturbation theory 
[26l |. Since the potential does not break 2D rotational 
symmetry, to calculate the correction to the ground state 
energy we only need to know the s-wave eigenvalues e„ 
and eigenfunctions -tjjn of the 2D harmonic oscillator: 



En = 2(2n+ 1), 



(A5) 



where L„ are the Laguerre polynomials. In calculating 
the matrix elements of the perturbation we use the iden- 
tity 



i=0 
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which follows by induction from the recurrence relation 
for the Laguerre polynomials [13]. After straightforward 
algebra we find 



2 Szcq 



(A7) 



Requiring the second term to be a small correction gives 
the condition Pr ^ 5.44; hence the expansion should be 
reliable up to Pr of order 1. Indeed, both the calculation 
of the coefficient of the 0^ term (~ 0.003) and comparison 
with numerics (see the end of this Appendix) show that 
Eq. (jA7p is still a good approximation even for /?,. > 1. 



2. Large /?, 

In the limit of large the kinetic term in Eq. ([7]) be- 
comes dominant, which physically correspond to an al- 
most free particle, and the (properly defined) potential 
U{p) can be treated as a shallow one. Since P{p) 1 
at large p, we define U{p) — P{p) — 1 and e = 1 — Cqi so 
that positive energies correspond to bound states. In two 
dimensions the ground state energy in a shallow poten- 
tial depends exponentially on the inverse mass [26|: e oc 
/Sr-e"'^^'', where the constant c cx | / cPpU{p)\~^ . This re- 
sult holds when the integral converges, while in our case 
the integral is logarithmically divergent, see Eq. (jA2[) . 
We can adapt to the present situation the derivation of 
the above estimate for e presented in Ref. [26|, §45. We 
show at the end of this appendix that the expression thus 
obtained agrees well with numerical calculations. 

We define 



y 



J/ 



and rewrite Eq. ([7]) as 



with 



U{v) = 1 + 
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Here we include the eigenvalue term (cx e in the original 
equation) into the potential energy U . In a small cir- 
cle around the origin of radius j/o the potential can be 
approximated as 



u{v) 



(All) 



In this region the wave function can be taken as approx- 
imately constant, ^2 = 1) and integrating both side of 
Eq. ([X9|) gives 



2/0*2(2/0) = -^?/o 



(A12) 



To estimate j/Qj we note that according to Eq. (jAip the 
length scale over which the potential varies appreciably 
near the origin is x '--^ I/cq, so we take ?/o = ys/Pr/cQ. 
For y 3> \fsJWr^ the potential V is approximately 



(A13) 



and Eq. (jA9p has as solution the modified Bessel function 
of imaginary order K^j ^^{y). We want to match its 
logarithmic derivative to the estimate in Eq. (|A12[) : to 
do so, we use the following approximate expression valid 
for y, 1/V^ < 1 m 



' I3r sin 



1 A y 

T= log - 



- IE 



(A14) 
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Direct inspection of the matching condition at yo suggests 
taking e in the form 



(A15) 



for some parameter b. Then the matching condition re- 
duces to 



1 



1 



and solving for h we find 



6=2 arctan 



(A16) 



(A17) 



Substituting this back into Eq. (jA15|) we finally arrive at 

£0 = 1- 4^e-^lB^2^p^^-2^^rct^n[cl2^] ^^^g) 



Note that, in agreement with Eq. (jA14|) . K^^^^{y) has 
an infinite number of zeros at (approximate) positions 
yn = 26^''"^^^'''^, n — 1,2,3,..., and no other zeros 
at y > yi. It is easy to check that yn < yo, so that the 
approximate wave function constructed in the course of 
this derivation has no zeros, as expected for the ground 
state. 



To check the accuracy of the estimates in Eqs. (jATj) 
and (jA18[) . we also solve Eq. ([7]) numerically. That is, we 
first find a numerical solution to Eq. (|4]) for f{x); then 
we find numerical estimates of the ground state energy of 
Eq. ([7]) for various Pr- Interpolation of these numerical 
results gives the solid curve in Fig. [21 where we also plot 
Eqs. (|T7l) and (jAlSp for comparison. Using Eq. (|T7| for 
/3r ^ 3 and Eq. (|A18p for Pr ^ ^ gives estimates that 
deviate less than 1% from the numerics - see the inset of 
Fig.m 
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